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We present the identification of a bright submillimetcr (submm) source, SMM J163555. 5+661300, 
detected in the lensing cluster Abell 2218, for which we have accurately determined the position using 
observations from the Submillimeter Array (SMA). The identified optical counterpart has a spectro- 
scopic redshift of z = 4.044 ± 0.001 if we attribute the single emission line detected at A = 6140A to 
Lyman-a. This redshift identification is in good agreement with the optical/near- infrared photometric 
redshift as well as the submm flux ratio Stso/Sgso ~ 1-6, the radio-submm flux ratio 6*1.4/5850 < 0.004, 
and the 24/xm to 850/xm flux ratio S^/Sgso < 0.005. Correcting for the gravitational lensing amplifi- 
cation of ~ 5.5, we find that the source has a far-infrared luminosity of 1.3 x 1O 12 L0, which implies 
a star formation rate of 230M Q yr~ 1 . This makes it the lowest-luminosity SMG known at z > 4 to 
date. Previous CO (4-3) emission line obserations yielded a non-detection, for which we derived an 
upper limit of the CO line luminosity of L co = 0.3 x 10 10 Kkm/spc~ 2 , which is not inconsistent with 
the L co — Lfir relation for starburst galaxies. The best fit model to the optical and near-infrared 
photometry give a stellar population with an age of 1.4 Gyr and a stellar mass of 1.6 x 10 10 M Q . The 
optical morphology is compact and in the source plane the galaxy has an extent of ~ 6 x 3 kpc with 
individual star forming knots of < 500 pc in size. J163556 is not resolved in the SMA data and we 
place a strict upper limit on the size of the starburst region of 8kpc x 3kpc, which implies a lower 
limit on the star formation rate surface density of 12 M Q yr _1 kpc 2 . The redshift of J163556 extends 
the redshift distribution of faint, lensed SMGs, and we find no evidence that these have a different 
redshift distribution than bright SMGs. 

Subject headings: galaxies: individual (SMM J163556+661300) — galaxies: high-redshift — submil- 
limeter 
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1. INTRODUCTION 

The redshift distribution of submillimeter galaxies 
(SMGs) is vital to our understanding of the nature of 
these galaxies and their role in general picture of galaxy 
formation and evolution. The first systematic identifica- 
tion and redshift determination of SMGs relied on the ra- 
dio counterp arts and was carried out by IChapman et al.l 
(|2003l 120051 ). While the radio and optical/near- infrared 
counterparts have much better determined positions 
compared to the original 15" resolution submillimeter ob- 
servations and thus is an efficient approach for identifying 
the underlying galaxy, the radio emission rapidly dims 
with the distance and thus is likely to introduce a bias 
against SMGs with redshift z > 3.5. Due to the negative 
if-correction (sub)mm observations offer a great oppor- 
tunity to search for galaxies beyond z > 4 and one would 
expect to find SMGs up to redsh ifts z ~ 6 — 7. Through 
modeling of the SMG population IChapman et all (2005) 
predict that about 10% of the population is at red- 
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shift z > 4. IWall et al.l (|2008f) suggested a bimodal 
redshift distribution, where the brightest SMGs domi- 
nate at the highest redshifts and the fainter SMGs dom- 
inate towards low er redshifts, as initially proposed by 
llvison et alj ((2002), equivalent to 'cosmic down-sizing'. 

To date, there has been no systematic samples 

of very high reds hift (z > 3 .5) S MGs, however 

iDannerbauer et all (|2002l I2004D and lYounger et all 
(|2007f ) have presented several high redshift candidates 
that were identified through (sub-)mm interferometry 
(with and without radio counterparts). Three SMGs 
with accurat e redshifts are COSMOS J100054+023436 
(z = 4.547; ICapak et al.l l2008t ISchinnerer etai] l2008h 
and GN20 and GN20.2 (both z = 4.05: iDaddi et all 
2008b). These three sources have CO redshifts provid- 
ing solid evidence for correct identification. The galaxy 
GOODS 850-5 (a.k.a. GN10) does not have an optical 
spectroscopic redshift, but based on photometric red- 
shifts two in dependent groups co nclude it is z > 4 , pos- 
sibly z ~ 6, (|Dannerbauer et alj|2008t IWang et al.ll2007l . 
I2009T ) . In a recent paper. IDaddi et al.l (|2009l ) presented 
a redshift based on a single millimeter wavelength emis- 
sion line, which they interpret as the CO(4-3) line plac- 
ing the galaxy at redshift z — 4.0424. Very recently, the 
highest redshift SMG known to date has been identified, 
namely the z = 4.76 galaxy L ESSJ033229.4-27561 9 in 
the LABOCA map of E-CDF-S (|Coppin et al.ll2009h . 

In this paper we present the identification of 
the optical counterpart of the submillimeter source 
SMMJ163555.5+661300 (henceforth J163556), which 
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was detected in th e deep SCUBA^ map of Abell 2218 
(jKneib et alJ 12004 iKnudsen et~afl l2006h with a lensed 
flux of £"850 = 11.3 mjy. It is important to note that this 
is a 10(7 SCUBA detection, rendering it a very solid de- 
tection and not a spurious source. Exploiting the strong 
gravitational lensing effect by the cluster potential, these 
SCUBA observations have yielded the largest number 
of faint submm galaxies detected in one field. Here we 
present multiwavelength followup and spectroscopic ob- 
servations. 

Throughout we will assume an fl — 0.3, A = 0.7 cos- 
mology with Hq — 70kms~ 1 Mpc -1 . With this cosmol- 
ogy, at z = 4 1" corresponds to 6.95 kpc. 

2. OBSERVATIONS 

2.1. Submillimeter interferometric observations 

SMM J 163556 was observed at the Submillimeter Array 
(|Ho et al.l|2004[ ) on Mauna Kea in compact configuration 
on UTC 05 May 2009 using 7 of the 8 SMA antennas. 
The resulting image created has uw-coverage from 13.7 
kA to 121.6 kA at the observed local oscillator frequency 
of 340 GHz. Observations were made in very dry con- 
ditions, with T225 ranging from 0.04 to 0.06. The full 4 
GHz (2 GHz in each sideband separated by 10 GHz) were 
combined to achieve a total rms noise of 1.9 mjy in the 
final map. 

The SMA data were calibrated using the MIR software 
package developed at Caltech and modified for the SMA. 
Gain calibration was performed using the nearby quasars 
1642+689, 1638+573, and 1849+670. Absolute flux cal- 
ibration was performed using real-time measurements of 
the system temperatures, with observations of Callisto, 
Ganymede, Uranus and Neptune to set the scale. Band- 
pass calibration was d one using Callisto . The data were 
imaged using Miriad (jSault et al.| [l995). where the re- 
sulting synthesized beam is 2.1" x 1.5" with a position 
angle P.A.= 39deg. To verify the astrometry, the posi- 
tion of 1638+573 was calibrated using 1642+689 as refer- 
ence sources; both calibrators have accurately known po- 
sitions known from the International Celestial Reference 
Frame. A fit to the Mi>-data of 1638+573 shows that it de- 
viates from the known radio position by AR.A. = 0.14" 
(0.10") and ADec. = 0.02" (0.005") in the lower (and up- 
per) sideband. As the distance between 1638+573 and 
1642+689 is 11.6 degrees, four times larger than that 
between 1642+689 and J163556, the fit provides a con- 
servative upper limit to the systematics uncertainties in 
the SMA astrometry. 

2.2. Radio observations 

Radio maps of A2218 at 1.4 GHz and 8.2 GHz were pre- 
sented bv lGarrett et alJ (|2005D . The depth of the 1.4 GHz 
data, which were obtained with the Westerbork Radio 
Synthesis Telescope, is la r.m.s. — 15/iJy, and the depth 
of the Very Large Array 8.2 GHz map is —6 juJy. J163556 
is not detected in the two maps, placing 3<7 upper limits 
of 5*1.4 < 45/iJy and Ss.2 < 18/iJy. 

2.3. Optical and near-infrared imaging 

7 Submillimetre Common-User Bolometer Array l|Holland et al.l 
119991) mounted at the James Clerk Maxwell Telescope (JCMT) 



High angular resolution space-based optical and near- 
infrared observations have been obtained with the Ad- 
vanced Camera for Surveys (ACS) and NICMOS on the 
Hubble Space Telescope for the programs 9292, 9452, 
9717, 10325, and 10504 with a total exposure time in 
the range from 1.6 to 5.5 hours. The data reduction 
of the se data is discussed in detail in Rich ard et al.l 
(2008). In summary the ACS data were reduced us - 
ing the multidrizzle software (Koe kemoer et all l2002'l. 
which also removes cosmic rays and bad pixels and com- 
bines the dithered frames to correct for camera distor- 
tions. The NICMOS data were reduced following the 
procedures given in the NICMOS data reduction hand- 
book, then bad pix els were flagged, an d subsequently us- 
ing the LACDSMIC ()van Dokkumll2001f ) IRAF procedures 
for cosmic ray rejection. 

Ground-based observations in the i-T-band were ob- 
tained using MOIRCS on Su baru (Multi-Object Inf raRed 
Camera and Spectrograph; llchikawa et alJ [2006;) as a 
part of t he large search of ext remely high redshift lensed 
galaxies (Ric hard et al.ll2~008f) . The integration time is 
—5 hr with good seeing conditions (0.3"-0.4") using 
dithered exposures of 50 seconds duration. The data 
were reduced using the MCSRED software package 8 an d 
the details on this is described in lRichard et a l. (2008). 

The astrometry of the optical and near-infrared images 
was obtained using the USNO and the 2MASS catalogs 
for the field and the dispersion gives an accuracy < 0.2". 

Optical and ne ar-infrared photometry w as performed 
using SExtractor (|Bertin fc Arnoutslll996l ) for a 0.5 arc- 
sec aperture for the HST data and with a standard aper- 
ture correction in the X-band. When undetected, we 
place 3cr upper limits. The photometry for the source 
z4, which is at a distance of 0.6 arcsec from the SCUBA 
position, is given in Ta ble [1] and all magnitudes are in the 
AB system (|0kelll974l ). The images are shown in Fig. [TJ 

2.4. Mid-infrared imaging 

The ultra-deep (10 hr/band) IRAC data of Abell 2218 
were obtained as part of the Spitzer GO-2 program PID 
20439 (PI: Eiichi Egami) while the deep (-3000s) MIPS 
24 /im data were taken as part of th e Spitzer GTO pro - 
gram PID 83 (PI: George Rieke). See lEgami et~aT1 (|2006l ) 
for the details of the IRAC and MIPS data reductions. 

J163556 is detected in all four IRAC bands (3.6, 4.5, 
5.8 and 8.0 /im) but not at 24 fim. A convolved HST im- 
age was used for determining the HST-IRAC colors. The 
convolved HST image was also used for the astrometry 
registration. The MIPS non-detection sets a 3cr upper 
limit of < 40 /*Jy at 24 /im. J163556 is not detected in 
the Inf rared Space Observato ry ISOCAM maps at 7 and 
15 ^m ([Metcalfe et al.ll2003D . We place 3cr upper limits 
of 63 and 150 /im. Jy, respectively. 

2.5. Optical spectroscopy 

We conducted deep multi-slit spectroscopy using the 
Low Resolution Imaging Spectrograph (LRIS; mounted 
at the 10m Ke ck Telescope on Mauna Kea, Hawaii; 
lOke et all Il995f ) of sources in rich galaxy cluster field 
A 2218. The data were obtained during the nights of 
June 30 and July 1, 2003, which had reasonable seeing, 

8 Available from http://www.naoj .org/staf f /ichi/MCSRED/mcsred.html 
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0.8", but were not fully photometric (with some cir- of continuum on the blue side. 



rus). We obtained a crude flux calibration of our ob- 
servations using Feige 67 and 110 as spectrophotometric 
standard stars. We observed the galaxy z4 (see Section 
[2~3]) for a total of 11.8 ks using the 400/8500 red grat- 
ing and a total of 10.6 ks using the 600/4000 blue grat- 
ings offering respectively a spectral dispersion of 0.63A 
pixel -1 in the blue and 1.86A pixel" 1 in the red. No flux 
is detected in the blue channel below 5600A, where the 
dichroic splitter is located. We reduced these data that 
was part of a multi-object mask using standard IRAF 
scripts. 

3. RESULTS 

We fitted a point source to the SMA data to determine 
the position of the J163556. When fitting to the uv-data 
we find point source at the position RA, Dec (J2000) 
= 16:35:55.67, 66:12:59.51 with positional errors of 0.15 
arcsec with a flux of 10.7±2.0mJy, and consistent results 
are found for the upper and lower sidebands, individually. 
When fitting to the synthesized map, we find RA, Dec 
(J2000) = 16:35:55.64±0.02, +66:12:59.44±0.14 and a 
flux of 10.7 ± 0.2 mjy, with similar results obtained using 
routines JMFIT from AIPS and imfit from MIRIAD. We 
note that this po sition is only an arc second from the 
SCUBA position (Kn udsen et al.ll2008l ). demonstrating 
that good positions can be obtained for sources in the 
coarse resolution maps when the signal-to-noise ratio is 
sufficiently high (in this case ~ 10). The flux determined 
by the SMA data agrees with the SCUBA flux, further 
verifying the source is unresolved in the SMA image. 

In Fig. [1] we show the SMA data overlaid on the optical 
data. The multiwavelength photometry for the underly- 
ing counterpart galaxy of J163556 is presented in TableHJ 
The optical/near- infrared position is determined in filter 
bands F625W to K. Because of the complex morphology 
of z4, we smoothed the ACS images with a 0.6" Gaus- 
sian kernel and the F110W band with a 0.5" Gaussian 
kernel before determining the position. The morphology 
varies between the different optical/near- infrared bands. 
We derive a mean optical/near-infrared position for z4 of 
RA,Dec(J2000) = 16:35:55.67±0.03, +66:12:59.42±0.08. 
The offset between this mean position and the SMA po- 
sition derived from the wt>-data is 0.09", and the uncer- 
tainty on the offset is 0.23" in R.A. and 0.17" in Dec, 
hence we conclude that optical and submm position agree 
well. 

Figure [2] shows the 2D and ID extraction of this object 
within the 6000-6700A window, where we have identified 
an (asymmetric) emission line centered at A = 6135A. 
Two independent reductions of the data were carried out 
producing similar results, and thus the emission line is 
most probably not spurious. A faint continuum is weakly 
visible red wards of the detected emission line. In light 
of other evidence discussed below, we believe this line to 
be Lyman-a, which correspond to a redshift z = 4.044 ± 
0.001, obtained through fitting to the emission line and 
a few absorption lines. We note these absorption lines 
are marginal and thus carry a low weight. A break in 
the continuum bluewards of the line as well as the lack 
of continuum emission in the blue part of the spectrum 
supports this identification. Alternatives at low redshift 
such as [OIII] at z = 0.23 or [Oil] at z = 0.65 seem less 
likely due to the absence of other lines and/or the lack 



4. DISCUSSION 
4.1. Redshift 

The ratio of the 450 (im and 850 /im fluxes of 1.6 
is strongly indicative of J163556 being a high redshift 
source with z > 3.5, assuming that the far-infrared spec- 
tral energy distribution (SED) is well-described by a 
modified black body (as illustrated in Fig. [3]). A high 
redshift is supported by the 1.4GHz to 850 /im flux ratio 
of < 0.004, where the depth of the radio data allows a 
constraint of z > 3. In Fig.[3]we show how the 1.4GHz to 
850/xm flux density ratio changes with redshift usi ng the 
relatio n fro mlCarilli fc Yuri (1 19991) . As discussed in lBlainl 
( 1999 1) and ICarilli fc Yun (1999), the radio-submm flux 
can be used as a redshift indicator, assuming that both 
the radio and submm emission is powered by star forma- 
tion. Additional evidence for J163556 being at high red- 
shift [z > 3.5) comes from th e 24/im to 850/im flux ratio 
5*24/^850 < 0.005 (e.g. Fig. 3 lPope et al.ll2006h . We note 
that the individual flux ratios in the radio, submm, and 
far-infrared, only provide crude redshift estimates and 
are ea ch known to have a large scatter (e.g. lCarilli fc Yunl 

l20onh . 

As the optical spectrum has only a single emission line, 
we ca lculate the photometri c redshift using the Hyperz 
code (Bolzonc lla et al.l l2000t ) and the photometry from 
the F435W (= B) band to the IRAC bands. We see an 
excess in the 3.6/im band (shown i n Fig.[H), which is most 
probably caused by Ha emission (|Yabe et alJ l2008). For 
redshifts z = 3.84 — 5.03 the Ha line will shift through 
the 3.6/im IRAC channel. We make a correction for the 
Ha emission deriving a line luminosity based on 10% of 
the star formation rate assuming a typical FIR luminos- 
ity for SMGs. Excluding the 3.6/^m data point, we find a 
photo- z of 4.67, while including the Ha corrected 3.6/im 
data point we find a photo- z of 4.77. This high photo- 
metric redshift strongly supports our identification of the 
emission line in the optical spectrum as Lya and not e.g. 
low redshift [Oil]. 

Adopting a redshift of 4, we estimate the gravitational 
lensin g magnification using t he cluster potential model 
from (lEHasdottir et all l2007f) as input to LENSTOOL 
(| Julio et alJl2007f ). We find a magnification factor of 
5.5 ± 1.9 9 and note that this does not change much for 
redshift z > 2. The spatial magnification is a factor 3.2 
in the direction tangential to the direction of the cluster 
center and 1.7 in the direction perpendicular. In Fig. [5] 
we show the reconstruction of the source plane image of 
J163556 using the WCS images. The morphology of the 
object is only slightly distorted, albeit all the scales are 
reduced by a factor of ~ 2. The different colors of the 
clumps are likely due to different reddening. 

The CO (4-3) line was not detected for J163556 in 
IRAM Plateau de Bure observations in the redshift range 
z = 4.035-4.082 (jKnudsen et al.ll2009h . Assuming a line 

9 We note that this value, although diff erent, is consistent wit h 
the magnification factor that we used in Knudsen ct al. (2009). 
Th e difference is a conseq uence of a recent revision of the model 
in ( Eh'asdottir et al. 2007), which aff e cts th is region of the cluster, 
and the fact that in I Knudsen et all ((2009) was determined for a 
single position, while here an extensive analysis to determine the 
uncertainty was carried out. Differences of 12%-20% are generally 
to be expected. 
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width of 400 km/s lKnudsen et all (|2009t ) placed an upper 
limit on the CO line luminosity of 0.3 x 10 10 Kkm/spc 2 
(using a lcnsing magnification of 5.5). This upper limit 
is not inconsistent with the L' co — Lfir relation for 
starburst galaxies 10 , hence the non-detection does not 
contradict the determined redshift z = 4.044. It is 
important to keep in mind that we do not know if 
the CO (4-3) transition is thermalized. While several 
bright SMGs have been found to have th ermalized CO 
at le ast up to the (3-2) transition (e.g. iTacconi et al.l 
2008), cases of Milky Way type molecular gas excitation 
conditions have been found for a z ~ 1.5 B zK galaxy 
(|Daddi et alj|2008at iDannerbauer et~aT1 l2~009) as well as 
subthcrmal excitation of the CO(2-l) and C O(5-4) in the 
submm bright extremely red objec t HR10 (jGreve et al.l 
[20031: iPapadopoulos fc Ivisonll200l . The CO (4-3) non- 
detection could also be explained by a velocity shift be- 
tween the Lyman-a and the CO larger than — 525km/s, 
which is the difference between z = 4.035 and 4.044 for 
CO (4-3). Shifts of ZL ya — zgo = — 500 km/s have been 
seen for other SMGs (Greve et al.ll2005l ). 

4.2. Spectral energy distribution 

Based on the 850/im and the 450/im data points, the 
far-infrared SED is well-described by a modified black- 
body spectrum, with a dust temperature ~ 35 K assum- 
ing /3 = 1.5, where /3 is f v cx v@ . The dust temperature 
is comparabl e to what is seen in other SMGs at lower 
redshift (e.g. iKovacs et ai]|2006l ). If the single optical 
emission line would not be Lyman-a, but [OIII] or [Oil] 
implying a low redshift, the implied dust temperature 
would be 9 — 11 K. In Fig. [6] we show the multiwave- 
length SED overlayed with a m odified blackbody spec- 
trum and the Arp220 SED from lSilva et~aTl (|1998h . The 
far-infrared luminosity derived for this redshift 4.044 is 
1-3^q - 3 x 10 12 L© after correcting for gravitational lens- 
ing and where the two uncertainties represent the up- 
per and lower limit from the magnification correction. 
J163556 is intrinsically a fainter object than the other 
z > 4 SM Gs, which have estimated FIR luminosities 
~ 10 13 L^ (iCapak et all 127)081 : ID addi et al.ll2008bL [20091 
iCoppin et al.l 120091 ^ Assuming a Salpeter initial mass 
function, solar abundance, a continuous burst for 10-100 
Myr, and that the dust reradiates all the bolometric lu- 
minosity, the far-infrared lumi nosity implies a star for- 
mation rate of ~ 230M Q yr- 1 (|Kennicuttlll998h . 

Assuming that the submm and radio emission from 
J163556 are correlated, we use the radio-su bmm f l ux ra - 
tio as function of redshift as discussed by iBlainl (|1999l ) 
to estimate the radio flux. We estimate Si. 4 ~ 40/iJy. 
Assuming that the radio SED is described by a power 
law, S v oc v a with a = —0.7, then S&.2 ~ 13/iJy. Both 
estimated fluxes are below the sensitivity levels of the 
WSRT and VLA maps. This implies that no radio-loud 
AGN is present and thus the submm and radio emis- 
sion is most likely powered by star formation rather than 
AGN activity. The absence of a dominating AGN or 
quasar is also supported by the l ack of an X-ray detec- 
tion in the Chandr a observations (jMachacek et alj|2002t 
iGovoni et "all 120041 ): in the 0.8 - 4keV band an upper 
limit of 1.24 x 10~ 14 erg/s/cm 2 is estimated for an area 

10 We note tha t the far-infrared luminosity given in 
IKnudsen et al.l 1 120091) unfortunately missed a factor (1 + z)" 1 . 



of 3 x 3 pixels. We note that this is the first z > 4 
SMG identification without radio counterpart, although 
we also point out that the radio observations are not as 
deep as those used for the identification of other z > 4 
SMGs. 

To determine the properties of J163556 we use hyperz 
and find that the optical and near-infrared photometry 
is best fit by a star formation history of an elliptical 
galaxy template with an age of ~ 1.4 Gyr and a red- 
dening of Ay = 0.8 mag. The rest-frame if- band ab- 
solute magnitude is -24.40. For this template the K- 
band luminosity would correspond to a stellar mass of 
8.7 x 10 10 M© without correcting for the gravitational 
lensing. Correcting for lensing magnification this cor- 
responds to ~ 1.6 x 10 10 M Q . The results of the fit, 
however, yield an age close to the age of the Universe 
at that redshift, ~ 1.5 Gyr. Additionally, we fix the red- 
shift to the spectroscopic redshift and use different star 
formation histories from a single burst to a continuous 
star formation. A single burst gives the best fit with a 
reduced \ 2 01 1-39, a reddening of Ay = 1.8 mag and 
very young age of 10-20 Myrs. Other star formation his- 
tories with an exponential decay or with continuous star 
formation yield a x 2 of 1.76, reddening of Ay = 2.0 mag, 
and a best fit age of 40 Myrs. 

We note that the Spitzer IRAC observations probe the 
light at rest-frame wavelengths of 0.71 — 1.6/im, which 
for an elliptical galaxy is dominated by the evolved stel- 
lar populations. However, for a starburst galaxy this 
could be dominated by young supergiants, as is seen in 
local ULIRGs. It is known that there is a degeneracy be- 
tween a mature stellar population with small or modest 
amount of extinction and a younger population with a 
much stronger extinction. A younger population would 
imply a lower stellar mass estimate and thus the best fit 
estimate can be considered an upper limit. 

Based on the F775W filter, we can compute the rest- 
frame UV continuum luminosity L(1500A) = 2.0 x 
10 22 W/Hz and derive an unlensed s tar formati o n rate 
(SFR) of - 28 M© yr" 1 using the iKennicuttl (Il998h 
empirical calibration (although certainly underestimated 
due to dust extinction). This result compared to the in- 
frared luminosity suggests that the vast majority of the 
young stars in this galaxy are obscured by dust and are 
undetectable in the rest-frame far-UV. 

4.3. Size and morphology 

The HST images give the best insight to the morphol- 
ogy of the galaxy. The resolution is the highest in the 
ACS F850LP (z'-band) image, which probes redshiftcd 
UV light (see Fig. Q}. The galaxy is seen as a very ir- 
regular system made of four knots. The total extent of 
the galaxy in the z'-band image is 1.8" x 1.1". Correct- 
ing for the lensing magnification, the source size is then 
0.9" x 0.4". For z — 4.044 this corresponds to a physical 
scale of ~ 6 kpc x 3 kpc. The largest star forming knot is 
unresolved in the ACS images and and we place an upper 
limit on the extent of about 0.07" or 500 pc in the source 
plane. In the if-band (which corresponds to rest-frame 
B-band) image the light originates primarily from the 
central part of the galaxy. We suggest that there already 
are old stars in place at the center of this galaxy. 

J163556 is detected as a point source in the SMA data. 
Correcting the beam shape for a gravitational lensing 
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amplification in the of 3.2 x 1.7 in the tangential and the 
radial direction, respectively, of the cluster potential well 
we find the beam to be 0.47" x 1.2" in the source plane. 
This is equivalent to the resolution that can be obtained 
with the SMA in extended configuration without gravi- 
tational lensing. Through a 2D-Gaussian fit to the SMA 
data we are not able to resolve the source concluding that 
it is point-like within the SMA beam. We place a strict 
upper limit on the size of the starburst region using the 
lensing corrected beam. In the z = 4.044 source plane 
this corresponds to 8kpc x 3kpc. For a star formation 
rate of ~ 230M Q yr _1 , this yields a lower limit for the 
star formation rate density of 12 M Q yr _1 kpc~ 2 . 

4.4. Extending the redshift distribution of faint SMGs 

The previously identified five z > 4 bright SMGs have 
far-infrared luminosities 5 — 10 times larger than J163556. 
For the six other reliably identified faint SMGs with 
■Sssoum < 2-mJy flFraver et all 120031: iKneib et "all 12004 
iBorys et al.ll2004t iKnudsen et alJl200a : Knudsen et al., 
in prep), the redshifts range from 1.03 to 2.9 with ma- 
jority of the sources found around z ~ 2.5. Due to the 
modest number of known faint SMGs, their redshift dis- 
tribution is poorly constrained and it has been proposed 
by several groups that the faint SMGs or low luminos- 
ity S MGs peak at low er redshifts than the bright SMGs 
(e.g. I Wall et al.ll200§ ). Our result for J163556 thus ex- 
tends the redshift distribution of faint SMGs to very high 
redshifts. 

In Fig. [7] we plot the redshift distribution of Ssso^m < 
2 mjy SMGs with spectroscopic redshifts in comparison 
to the reds hift distribution for ra dio identified, blank 
field SMGs (|Chapman et al.ll2005T) . Between redshift 1 
and 8, due to the negative fc-correction in the submm, 
the observed flux density stays essentially constant for a 
given far-infrared luminosity. As several SMGs have their 
far-infrared luminosity determined from single submm 
fluxes, we use the flux density rather than the luminosity 
for our simple sample comparison. As the radio identified 
SMGs are expected to represent the general SMG popu- 
lation, but with a possible bias against reds hifts z > 3.5, 
we also include the five z > 4 bright SMGs (ICapak et al.1 
[2008tlCopt)in et al.ll200l iDaddi et al.ll2008E I2009D . Us- 
ing the non-parametric ra nk sum test known as the 
Mann- Whitney U test (e.g. iBabu fc F eigclson ll996l ). we 
find that there is only a 1.3a (1.6<r) level of significance 
that the faint SMGs have a different redshift distribu- 
tion than the Chapman et al. sample + the z > 4 SMGs 
(the Chapman et al. sample alone) . Removing the z < 1 
tail of the Chapman et al. sample, as the negative k- 
correction does not apply in that redshift range, we find 
a consistent result. In other words, based on this small 
number of faint SMGs we find that it is possible that the 
less luminous SMGs have a similar redshift distribution 
as the more luminous SMGs. 

With the model pr edictions for the red s hift distri- 
butio n of SMGs (e.g. IChapman et al l 120051: IWall et all 
2008f) in most submm surveys (e.g., ISmail et al.1 120021: 
Webb et all 120031 : ICoppin et all 120061 : IKnudsen et all 
200l " at least a few sources are expected to have a 
z > 4. As z ~ 2 massive, old evolved galaxies (> 
10 11 M Q , > 2Gyr) are observed with a space density o f 
- 10- 4 Mp C - 3 Ce.g. lDaddi etal]l2005t IRong et al.ll200l . 
These galaxies probably had their dominant stellar pop- 
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ulation formed in a single starburst at r edshift z > 4 (e.g. 
iDaddi et alJl2005HStockton et al.ll200l . The most likely 
progenitor of these are the massive starbursts that we ob- 
serve in SMGs. If the large SFRs of - 1000 M yr" 1 can 
be sustained over a period of ~ 100 Myr it is possible t o 
build such large stellar masses (e.g. ICoppin et afl [2009). 
Thus identification of > 4 SMGs is vital for our under- 
standing of z ~ 2 massive, old galaxies. 

To reliably identify z > 4 SMGs, deep data at all 
wavelengths is necessary. Though the strongest, initial 
probe comes from the ratio between 450 /im and 850 /im 
fluxes, or possibly 350 /im and 850 /im, where this ra- 
tio is expected to be < 2. In most cases this requires 
very deep 450 /im or 350 /im observations with a la rms 
< 5mJy/beam. Thus the possibilities for systematic 
searches for very high redshift SMGs rely on the suc- 
cess of the future high-sensitive submillimeter cameras 
such as SCUBA-2 and space telescopes such as Herschel 
to provide such observations. 

5. CONCLUSIONS 

We report the identification of the gravitationally 
lensed SMG SMMJ163555.5+661300 (J163556), which 
has a lensing corrected far- infrared luminosity of 1.3 x 
10 12 L Q . Using SMA interferometric observations, we 
have identified the underlying optical counterpart. The 
SMA flux agrees well with the SCUBA flux of 11.3 mjy. 

• From optical spectroscopy from the Keck Tele- 
scope, we detect an asymmetric emission line at 
6135 A and no emission on the blue side. We iden- 
tify this emission line as the Lymana line, which 
implies a redshift of 4.0441 ±0.001. Using pho- 
tometric redshift calculations we rule out that the 
single emission line could be [Oil] or [OIII], which 
would otherwise have implied a redshift of 0.25 and 
0.61, respectively. This is also strongly supported 
by the flux ratios at submm and radio wavelengths. 

• J163556 is not resolved in the SMA data, and we 
place a strict upper limit on the starburst size of 
8kpc x 3kpc. From the derived SFR of 230 M Q yr" 1 
we estimate a lower limit on the SFR surface den- 
sity of 12M Q yr _1 kpc~ 2 . 

• With an unlensed flux of ~ 2mJy, J163556 is 
fainter than the classical (bright) SMGs that are 
often classified at Ssso^m > 5mJy. It is thus so 
far the most distant known faint SMG. This result 
extents the redshift distribution of the faint SMGs 
that are typically only found in gravitational lens- 
ing surveys to look-back times equivalent to the 
bright SMGs. 

• The derived upper limit for the CO (4-3) line lumi- 
nosity of 0.3 x 10 10 K km/s pc~ 2 is not inconsistent 
with the L C q — Lfir relation for starburst galax- 
ies. The CO(4-3) non-detection is possibly caused 
by subthermal excitation or a large velocity shift. 

• Based on the small number of reliably identi- 
fied faint, lensed SMGs with spectroscopic red- 
shift, we do not find evidence that they have a 
different redshift distribution than that of bright 
SMGs as traced by the radio identified SMGs from 
IChapman etaH (|2005f) . 
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Fig. 1. — Top panel: The F850LP (z') image overlayed with the SMA contours (white; in steps of lu rms. starting at 2a). The black 
circle shows the SCUBA beam (FWHM of 15"). The black rectangle shows the position of the slit used for optical spectroscopy. Middle 
panel: Mosaic of the optical and near-infrared images of J163556. The size of each image is ~4"x4". Bottom panel: Is a ~ 2" X 2" 
true-color image of the F555, F775, and F850LP overlayed with left: X-band contours and right: SMA contours, where the plus sign 
indicates the submm position as measured with the SMA (the plus is 0.2" X 0.2"). 
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Fig. 2. — Optical spectroscopic identification of J163556. (Top) The 2D sky subtracted spectrum that shows a faint emission line. It has 
been smoothed using a Gaussian kernel to visually enhance the emission. (Bottom) The ID sky subtracted spectrum. The hashed vertical 
lines indicate the location of bright sky lines. The best fit to the Lya emission as well as several (marginal) absorption lines yield a redshift 
of z = 4.044 ± 0.001 ; the dotted vertical lines show the position of the Lyo line and the absorption lines. The right panel shows a zoom in 
on the Lyct line. 




Fig. 3. — The predicted 450^tm-850/xm ratio and 1.4GHz-850^tm flux ratio as function of redshift. The former is c alculated assuming a 
modified blackbody SED (with T going from 30 to 50K from left to right) and the latter is based on the relation from Carilli & Yun (1999) 
(with a = —0.8 as the left curve). The horizontal dotted line shows the flux ratio and flux ratio limit, respectively. The vertical solid line 
indicates z = 4.044, while the vertical dashed lines show the redshift 0.23 and 0.65. 
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Fig. 4. — The SED from the F435W to 8.0/^m photometry overlayed with the best-fit spectrum from Hyperz from the photo- z modeling of 
the optical counterpart of J163556. The blue 3.6^tm shows the photometry as measured, and the black 3.6/^m data point has been corrected 
for He* emission line contribution. The blue spectrum shows the fit when including the Ha corrected 3.6/xm data point {z p ^ = 4.77) and 
the red spectrum shows the fit when excluding the 3.6/^m photometry from the photo-z estimate = 4.67). The insert shows the output 
probability distribution for the photometric redshift from Hyperz, where the solid line corresponds to the fit including the corrected 3.6/jm 
photometry and the dashed line corresponds to the fit with 3.6/^m excluded. 
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Fig. 5. — Reconstruction of the source plane image of J163556 shown using WCS images in the filters 625W, 775W, and 850LP. 
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Fig. 6. — The SED of J163556, overlayed with a modified blackb ody with = 30, 35,40 K (dashed, solid, dotted, respectively) and 
j3 = 1.5, as well as the SED model for Arp220 from Silva et al. (1998) scaled to the submm points of J163556. 








Fig. 7. — The redshift distribution for faint and bright SMGs. The solid (purple) histogram shows the redshift distribution of 
the reliably identified lensed SMGs with Ssso^tm < 2 mjy. The long dashed (green) histogram sho ws the redsh i ft dis t ribution from 
IChapma n et al ] II2005P), and the d otted (red) histogram shows the same including the five redshifts from Capak et al. (2008); Coppin et al. 
H2009l ): IDaddi et al.l (|2008bl . I2009T ) . 
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TABLE 1 
Photometry for J163556. 
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